PROCEEDINGS 

OF SCIENCE 



Predictions from Lattice QCD 



FERMILAB-CONF-05-428-T 



in 
o 
o 

(N 

Or 
<D 

o 

m 

(N 
> 

o 
in 
o 



X 



A.S. Kronfeld* fl , I.F. Allison 6 , C. Aubin cJ , C. Bernard^, C.T.H. Davies 6 , C. DeTar c , 
M. Di Pierrot E.D. Freeland#, Steven Gottlieb' 1 , A. Gray', E. Gregory^, U.M. Helled, 
J.E. Hetrick', A.X. El-Khadra ", L. Levkova /! , P.B. Mackenzie ', F. Maresca e , 
D. Menscher"', M. Nobes" , M. Okamoto M.B. Oktay m , J. Osborn e , D. Rennet, 
J.N. Simone ', R. Sugar'', D. Toussaint', H.D. Trottier 

E-mail: ask@ f nal . gov] 



"Fermi National Accelerator Laboratory, Batavia, Illinois 60510, USA 

b Department of Physics and Astronomy, Glasgow University, Glasgow, Scotland G12 8QQ, UK 

c Physics Department, Columbia University, New York, New York 10027, USA 

d Department of Physics, Washington University, St. Louis, Missouri 63130, USA 

e Physics Department, University of Utah, Salt Lake City, Utah 841 12, USA 

f School of Computer Science, Telecommunications and Information Systems, DePaul University, 

Chicago, Illinois 60604, USA 

g Liberal Arts Department, School of the Art Institute, Chicago, Illinois 60603, USA 
^Department of Physics, Indiana University, Bloomington, Indiana 47405, USA 
'Department of Physics, The Ohio State University, Columbus, Ohio 43210, USA 
1 Department of Physics, University of Arizona, Tucson, Arizona 85721, USA 
American Physical Society, Ridge, New York 11961, USA 

1 Physics Department, University of the Pacific, Stockton, California 95211, USA 

'"Physics Department, University of Illinois, Urbana, Illinois 61801, USA 

"Laboratory of Elementary-Particle Physics, Cornell University, Ithaca, New York 14853, USA 

"Physics Department, Simon Fraser University, Burnaby, British Columbia V5A 1S6, Canada 

p Department of Physics, University of California, Santa Barbara, California 93106, USA 

Fermilab Lattice, MILC, and HPQCD Collaborations 

In the past year, we calculated with lattice QCD three quantities that were unknown or poorly 
known. They are the q 2 dependence of the form factor in semileptonic D — » Klv decay, the 
decay constant of the D meson, and the mass of the B c meson. In this talk, we summarize these 
calculations, with emphasis on their (subsequent) confirmation by experiments. 
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1. Introduction and Background 

In recent years, lattice QCD has reached the stage where many calculations of hadron masses, 
mass splittings, and operator matrix elements agree with experimental measurements. The key 
has been the inclusion of sea quarks. The progress has been especially striking [[Tj] when the light 
quarks (sea and valence) are implemented as staggered quarks, with an improved action. 

Some of the ingredients of these calculations are controversial. Staggered quarks come in four 
tastes, three of which must be removed to obtain each individual flavor. For sea quarks, this is 
done by taking the fourth root of the fermion determinant; for valence quarks, by projecting onto 
the desired taste sector. Furthermore chiral perturbation theory must be modified [^j. Although 
evidence for the validity of these "tricks" is slowly accumulating, a proof remains at large |j3p. In 
addition, much of the success of Ref. \^ comes from hadrons with heavy quarks. Although debate 
on heavy quarks in lattice QCD seems to have subsided, checks are still useful. 

In this paper, we discuss three calculations, with emphasis on their subsequent experimental 
confirmation. They are the normalization and g 2 -dependence of the D — > Klv form factor; the 
decay constants of the D + and D s mesons; and the mass of the B c meson. Each tests a somewhat 
different combination of the ingredients, and the following table gives an informal guide: 



calculation 


light sea 


light valence 


heavy 


semilep tonic f+(q 2 ) 


** 


** 


** 


leptonic fo 


** 


*** 


** 


B c mass 


** 




*** 



The chiral extrapolation, which is more sensitive to valence quarks than sea quarks, turned out to 
be more important for the decay constant than the form factor. The B c meson has no light va- 
lence quarks at all, but one should expect an accurate calculation only if heavy-quark discretization 
effects are under control. 

Successful predictions are, of course, not a substitute for a proof. They are still useful. Even 
if the experts are confident of all the elements of their numerical calculations, non-experts are 
interested in an end-to-end check [Qj. The quantities discussed here are ideal candidates: they are 
straightforward to compute; the first "good" experimental measurements were not expected until 
this year; and new physics is unlikely to contribute significantly. 

2. Semileptonic D Decays 

Semileptonic decays such as D — > Klv are mediated by electroweak vector currents. The 
matrix element (K\V^\D) is parametrized by form factors. For a vector current there are two, but 
experimentally only the one called f+(q 2 ) is accessible; the rate from the other one, fo(q 2 ), is 
suppressed by ntj. Here q 2 is the momentum transferred to the lepton-neutrino system, falling in 
the range < q 2 < q 2 msx = (mj) — mx) 2 - In lattice QCD, discretization effects are smallest when 
the momentum p of the kaon is small, and then q 2 is not too far from q^^. 

Experiments usually measure the branching fraction and quote the normalization /+(0), after 
making assumptions about the q 2 dependence. While our results were still preliminary [Q], exper- 
imental results came out for the normalization of D — > Klv [|]] and D — > %lv [j7|]. The agreement 
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with our final results [g| is excellent. For example, we find /+ *(0) = 0.73(3) (7) [§] while BES 
measures f+~* K (0) = 0.78(5) [gj. Our calculations of the normalization are also consistent with 
the soft pion theorem, which states /o(<?max) = Id/ fx- 

In principle, the shape of the form factors can be computed directly in lattice QCD. In practice, 
we calculated at a few values of p and used the Becirevic-Kaidalov (BK) form []|] to fix the full q 2 
dependence of / + and /q. Then the normalization of / + comes mainly from /o through a kinematic 
constraint /+(0) = /o(0). The BK Ansatz and calculations near <7^ ax determine the shape. It was 
important, therefore, to measure the q 2 dependence experimentally. In photoproduction of charm 
off fixed nuclear targets, the FOCUS Collaboration was able to collect high enough statistics to 
trace out the q 2 distribution of the decay Jl0|]. This setup does not yield an absolutely normalized 
branching ratio, so one is left to compare f+{q 2 )/ '/+(0). 

In Fig. |l| we plot our result for f+(q 2 )/ /+(0) vs. q 2 /m|>». The errors from /+(0) must be 
propagated to non-zero q 2 , so for / + (q 2 ) // + (0) the errors grow with q 2 . Figure [j] shows 1 - o bands 
of statistical (orange) and all uncertainties (yellow) added in quadrature [JXTJ] . As one can see, the 
q 2 dependence of lattice QCD (curve and error band) and experiment (points) agree excellently, 
although the uncertainties are still several per cent. 



D — > KN 
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lo (statistical) 
la (stat + syst) 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 
2 2 
q lm D , 



Figure 1: Shape of form factor f+(q 2 )/ /+(0) vs. q 2 /m^*, compared with experiment [[To|], 



3. Leptonic D Decays 

We also computed the hadronic matrix element for the leptonic decay of charmed mesons, f D + 
and fo s . The first (experimental) measurements of f D + appeared in 2004, with three events from 
BES @ and eight from CLEO [jT3|]. Neither provides a stringent test of QCD, but CLEO-c was 



just starting its run and promised 5-8 times higher statistics by the Summer 2005 Lepton-Photon 



Symposium [gp. At Lattice 2004 [|14J, we presented preliminary results for f D +, based on one 
lattice spacing, a « 0.125 fm. Our aim was to extend the running to two other lattice spacings 
and, of course, to improve our understanding of other aspects of the calculation, such as the chiral 



extrapolation. Details are given in the ensuing publication [15]. We find 

f D + = 201±3±17MeV, (3.1) 

where the first error is from finite Monte Carlo statistics, the second is a sum in quadrature of 
several systematics. A conservative (but not naive) estimate of heavy-quark discretizations effects, 
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as discussed in Ref. [jig], is the second largest (largest) systematic on f D + C/b s ). A few days after 



our paper was posted on the arXiv, CLEO-c announced its new measurement Q17[ ] 

f D+ = 223 ±17 ±3 MeV, 



(3.2) 



based on 47 ± 8 events. At the 1-a level, the agreement between Eqs. pA\ > and ( fT2j ) is fine. One 
should keep in mind that the experiment actually determines |V c ^|/o+. CLEO-c [|17|] assumes that 
\Vcd \ = \Vus\ and uses a recent average of \V US \ from semileptonic K decay. 

It is interesting to look at the rif dependence of fj) s , shown in Fig. ||(a). Of course, quenched 
results vary widely, but we show one [ p~8[ ] earned out with similar choices for heavy quarks, renor- 
malization factors, etc. One sees a trend of fo s to increase with n/. A similar comparison of f D +, 
in Fig. ||(b), is less instructive, because the chiral extrapolations in Refs. [18, 19] stalled at large 
quark masses and are, hence, less reliable than in the present work. 
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Figure 2: Dependence of (a) fo, and (b) f D + on the number tif of sea flavors. Quenched (rif = 0) [|l8|]; 
rif = 2 jl9|]; rif = 3 jl^], Solid (dashed) error bars are statistical (statistical+systematic). 



4. Mass of the B c Meson 

The pseudoscalar meson is the lowest-lying bound state of a charmed quark and a b quark. 
CDF [ pQ| ] first observed it during Run I of the Tevatron in the semileptonic decay B+ — > J /yl + v. 
During Run II, D0 has confirmed the discovery in the same mode [Mj. Because the neutrino is 
undetected, the mass resolution in semileptonic modes is poor, ±(300-400) MeV. Now, however, 
the upgraded detectors are able to reconstruct hadronic modes, such as S+ — > J/yn + , which give 



much much better precision on mg c [22]. 



At Lattice 2004 we presented results in nearly final form [23], and posted the final results on 
the arXiv in mid-November [ |24| ] : 

m Bc = 6304 ± 12+^ MeV, (4.1) 
where the last error is a rough estimate of residual heavy-quark discretization effects. Soon after- 



wards, CDF announced a precise mass measurement. They find [25] 



m Bc = 6287 ± 5 MeV, (4.2) 

which agrees with Eq. ([O]) at slightly more than 1-a. 

Two comments are in order. First, the agreement at the gross level of the calculation with 



experiment shows that discretization effects are well under control with lattice NRQCD [27] and the 



206/4 



Predictions f mm Lattice QCD 



A.S. Kronfeld 



Fermilab method [28]. Of course, this follows from the careful application of effective field theories 
for heavy quarks [|^, 30]. Indeed, as seen in Fig. ||(a), almost no lattice spacing dependence is seen 
in the splitting A^x = «b c — {fh^ + rax) /2 that is at the crux of the calculation [26]. Moreover, it is 
striking how much the splitting A^x changes when sea quarks are included. Figure ||(b) compares 



Eq. ( J4. 1[ ) with an old quenched calculation [|2g] (and the measurement Q25Q). The solid error bar 
shows the non-quenching errors, and the dashed includes the estimate of the quenching error. The 
inclusion of sea quarks has reduced the splitting by a factor of three or four, bringing an essentially 
discrepant result into agreement. 
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Figure 3: (a) Dependence of the splitting A^x on the lattice spacing a. (b) Comparison of the quenched [gf 
n f = 2 + 1 [^4|, and experimental p5j ] values of iiib c \ the dashed line denotes the baseline {illy + m^)/2. 



5. Conclusions 



In the past year, three lattice-QCD calculations have been confirmed by experiment. FOCUS 
[10] confirmed the -dependence of the D — > Klv form factor CLEO-c [17] confirmed the D- 
meson decay constant [JT3J] ; and CDF [ p5| ] confirmed the mass of the B c meson [24]. To obtain these 
results it is essential to have heavy-quark discretization effects under control, as one expects from 
theoretical foundations [27, 28, 29, Furthermore, the comparison of quenched QCD, QCD 



with 2+1 staggered flavors, and experiment shows that sea quarks are needed to obtain agreement, 
and that staggered quarks (in these cases) capture the needed effect. 

This work has been supported in part by the U.S. National Science Foundation, the Office of 
Science of the U.S. Department of Energy (DOE), and the U.K. Particle Physics and Astronomy 
Research Council. Fermilab is operated by Universities Research Association Inc., under contract 
with the DOE. 
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